Unified approach to structure factors and neutrino processes in nucleon matter 
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We present a unified approach to neutrino processes in nucleon matter based on Landau's theory 
of Fermi liquids that includes one- and two-quasiparticle-quasihole pair states as well as mean-field 
effects. We show how rates of neutrino processes involving two nucleons may be calculated in terms 
of the collision integral in the Landau transport equation for quasiparticles. Using a relaxation time 
approximation, we solve the transport equation for density and spin-density fluctuations and derive 
a general form for the response functions. We apply our approach to neutral-current processes in 
neutron matter, where the spin response function is crucial for calculations of neutrino elastic and in- 
elastic scattering, neutrino-pair bremsstrahlung and absorption from strongly-interacting nucleons. 
We calculate the relaxation rates using modern nuclear interactions and including many-body con- 
tributions, and find that rates of neutrino processes are reduced compared with estimates based on 
the one-pion exchange interaction, which is used in current simulations of core-collapse supernovae. 

PACS numbers: 97.60.Bw, 26.50.-fx, 95.30.Cq, 26.60.-c 



I. INTRODUCTION 



Neutrino emission, absorption and scattering processes in nucleon matter play a crucial role for the physics of 
stellar collapse, supernova explosions and neutron stars [l], Q • Since the leptons in these processes interact weakly, the 
neutrino rates can be expressed compactly in terms of the response of nuclear matter to axial and vector probes. In 
many situations, the axial response is the more important, and in this paper we concentrate on this case, which for a 
system of nonrelativistic nucleons amounts to the spin or spin-isospin response. These responses have been calculated 
' by a number of groups [!, 0, H, 0] allowing for single nucleon quasiparticle-quasihole pair states. 1 However, this is 
insufficient for rates of neutrino processes involving two nucleons, such as neutrino-pair bremsstrahlung and absorption, 
and modified Urea reactions, in which two particle-hole pair states are necessary. The possible importance of two 
particle-hole pair states for neutrino inelastic scattering,in_particular for energy exchange and the formation of the 
neutrino spectra, has been emphasized by Raffelt et al. @, [3, LLC} • Bounds on the magnitude of the two particle- hole 
pair weight have been investigated in Ref. [ll[ and it has been shown how the two-pair response is directly related to 
the collision term in Landau's transport equation for quasiparticles [l2j |. 

Noncentral contributions to nuclear interactions, such as tensor forces from pion exchanges and spin-orbit forces, 
are essential for the two particle-hole pair response, as is clear from calculations of neutrino-pair bremsstrahlung 
and the modified Urea processes [l3j and from general considerations based on conservation laws [ll| . Neutrino-pair 
bremsstrahlung and absorption change the number of neutrinos and are key for equilibrating muon and tau neutrino 
number densities in supernovae. The standard rates for bremsstrahlung are based on the one-pion exchange model 
for nucleon-nucleon interactions [l3[ (in the context of supernovae, see for example Ref. [lj3]). This is a reasonable 
starting point, since it represents the long-range part and the leading noncentral contribution in chiral effective field 
theory for nuclear forces [lj|- However, the tensor force from pion exchange is singular at short distances, which in free 
space requires iteration in the spin-triplet channels [l5j . In addition, subleading noncentral contributions to nuclear 
interactions are important for reproducing nucleon-nucleon scattering for the relevant channels and energies [l6l |. 

The aim of this paper is to give a unified treatment of neutrino processes that includes one- and two-particle-hole 
pair states as well as mean-field (Fermi liquid) effects consistently, and to present improved rate calculations of these 
processes based on modern nuclear interactions beyond one-pion exchange and including many-body contributions. A 
convenient framework for doing this is Landau's theory of normal Fermi liquids. This work represents an extension of 
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1 The basic single-particle-like excitations we work with are quasiparticles and quasiholes that have properties quantitatively different 
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Ref. |12j, which included two particle- hole pair states only in leading order using diagrammatic perturbation theory. 
Here wc shall use the quasiparticle transport equation. This provides a useful framework for understanding the 
basic physics and for making detailed calculations. In this paper, we focus on neutral-current processes in normal 
(nonsupcrfluid) neutron matter. We leave for future work the application to mixtures of neutron and protons, charged- 
current reactions, and the extension to supcrfluid phases. 

This paper is organized as follows. Section [ill gives an introduction to neutrino processes and the dynamical 
structure factors. In Sect. IIII1 we discuss Landau Fermi- liquid theory, show that it represents a useful effective 
theory for neutrino processes in nucleon matter, and introduce the transport equation for quasiparticles. Using a 
relaxation time approximation, we solve the transport equation for density and spin-density fluctuations and derive a 
general form for the response functions in Sect. IIVI The response function includes contributions from one-particle- 
hole pair (corresponding to elastic scattering of neutrinos from nucleons) and two-particle-hole pair states (which 
enter calculations of inelastic scattering, and neutrino-pair bremsstrahlung and absorption). In Sect.fVl we calculate 
the appropriate relaxation times for the one-pion exchange interaction and for a general operator representation of 
the quasiparticle scattering amplitude. We present results in Sect. I VII based on modern nuclear interactions and 
including many-body contributions, and contrast these with rates obtained using the one-pion exchange interaction, 
which is typically used in supernova simulations. Finally, we assess the significance of the improved treatment of 
nuclear interactions for neutrino mean free paths, energy loss and energy transfer in supernovae. We summarize the 
improvements and conclude in Sect. IVIIl 



II. NEUTRINO PROCESSES AND DYNAMICAL STRUCTURE FACTORS 



For neutral-current processes, the weak interaction Lagrangian density for low-energy probes is given by 

c{x)=^{x)nx), (i) 

where Gf = 1.166 x 10~ 5 GeV~ 2 is the Fermi coupling constant and the weak neutral currents are l^(x) for leptons 
and j/j,(x) for hadrons. The neutrino contribution to the leptonic current is 

= -7s)i>v , (2) 

and for nonrelativistic nucleons the hadronic current can be written as 

j^(x) = ?/' n 7m( c V - Ca75)^n ~ 0n(Cv<V) _ CaV' (Ji)<pN , (3) 

where ip v are neutrino fields, V'n nucleon Dirac fields, (f>N nonrelativistic nucleon spinors, and Ui Pauli matrices. The 
neutral-current vector coupling constant is Cy = — 1/2 for neutrons and Cy = 1/2 — 2 sin 2 9w ~ for protons, Ca is 
the axial-vector coupling, Ca = — ffa/2 = —1.26/2 for neutrons and Ca = Sa/2 for protons. While the vector current 
is conserved, the axial coupling can be modified in a many-body system. As a result, one may expect a reduction of 
g a for a nucleon quasiparticle by 5 — 10% in neutron matter and 10 — 20% in symmetric nuclear matter [l7l [l8j . 

Consider neutrinos with incoming energy lo\ and momentum qi that scatter from nuclear matter to a final state 
with energy u>2 and momentum q2 . Since neutrinos interact weakly, the rate for neutrino scattering can be expressed 
in terms of the dynamical structure factors for vector and axial responses of the nuclear medium [l|, Because 
neutron velocities in neutron matter at the densities of interest are nonrelativistic, these reduce to the density and 
spin responses, which are decoupled if the system is not magnetically polarized. 

The dynamical structure factors depend on the energy and momentum transferred to the system, lo = u)% — u>2 and 
q = qi — q2 , and are defined for the density response by [J, [l!| 

11 1 Z" 00 

S v (o;,q) = — z - 7F Imx(w,q) = - dt e*«" (n(t, q) n(0, -q)> , (4) 

and for the spin response by 

11 If 00 
S A>ij (w,q) = — , T ImXjj^q) = - / dte lut ( Si (t, q) Sj -(0, -q)> , (5) 

where n denotes the neutron number density, T is the temperature, s = <f>*(j<f> the spin density, and x(w,q) and 
Xij(w,q) are the density-density and spin-density-spin-density response functions, respectively. We use units with 

h = c = fcs = l. 
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In the long-wavelength limit, q — > 0, the spin response is in the direction of the applied magnetic field, hence Xij — 
for i ^ j. This is not the case at non-zero q, and the transverse and longitudinal spin responses differ However, 
for neutrino processes in supernovae and neutron stars, the momentum transfers are small compared with typical 
momenta of the nucleons, such as the Fermi momentum or the inverse Compton wavelength, and therefore the spin 
response is essentially diagonal, 

Xtj~XtrSij and SA,ij « Sa Sij ■ (6) 

The transition probability T(Qi, Q2) for a neutrino with energy and momentum Qi = {oj\, qi) to scatter to a state 
Qi = (t^2,q2) is fully determined by the density and spin response functions (see for example Refs. [11,131), 



T(Q 1 ,Q 2 ) = 2nnGl 



Cl (1 + cos 9) S v (u>, q) + C\ (3 - cos 9) S A (fi>, q) 



(7) 



where 9 = arccos(qi • q 2 ) is the scattering angle. The rate for bremsstrahlung of a neutrino with four- momentum Q2 
and an antineutrino with four-momentum Qi is given by T(—Q\, Q2), and for absorption of a neutrino with Qi and 
antineutrino with Q9 by T(Qi, —Q?)- We note that Eq. ([7]) neglects corrections of order u>/m from weak magnetism 
and other effects [20(. In terms of the transition probability, the rate of change of the neutrino occupation number 
n qi due to interaction with the nuclear medium is given by 



driq 1 f dq 



dt J (2tt 



r(Q2 ! Qi)nq 2 (l-nq 1 )-r(Qi,Q 2 )n qi (l-n q2 )+r(-Q2,Qi)(l-nq 1 )(l-«q 2 )-r(Qi,-Q 2 )nq 1 riq 2 



. ( 8 ) 

where n qi is the antineutrino occupation number. The four terms correspond to "in-scattering" , "out-scattering" , 
neutrino-pair bremsstrahlung and absorption, respectively. These differ only by the kinematics in the dynamical 
structure factors. 

III. LANDAU FERMI-LIQUID THEORY AND QUASIPARTICLE TRANSPORT EQUATION 

In supernovae and neutron stars, the neutrino energies are typically uii , u>2 < 30 MeV. The corresponding neutrino 
momenta <?i,<72 J$ 0.15 fm _1 are therefore small compared with the momenta of neutrons, which are of the order 
of the Fermi momentum fcp ~ 1.0 fm -1 for densities n ~ no/10. Here, no = 0.16fm -3 or po — 2.8 x 10 14 gcm~ 3 
denotes the saturation density of symmetric nuclear matter. Consequently, it is a good first approximation to work 
only to lowest order in the neutrino momenta. In addition, we focus on situations when the temperature is small 
compared with the Fermi energy of neutrons. This is the regime in which Landau's theory of normal Fermi liquids 
may be used [2ll |22|. Landau theory provides a clear separation between long- wavelength, low- frequency degrees of 
freedom, which are treated explicitly, and short-wavelength, high-frequency ones, whose effects are included in low- 
energy constants that incorporate the renormalization of matrix elements of currents and interparticle interactions. 
Another strength of Landau Fermi-liquid theory is that it brings out clearly the role played by conservation laws. 
Low-temperature expansions for Fermi liquids are often useful for T/e-p = I/77 < 1/tt. We therefore expect our results 
to be reasonable for degeneracy parameters ry > 3, which is typically valid for the relevant densities in supernovae and 
neutron stars. 

Nucleon matter differs from liquid 3 He, the prototype Fermi liquid, in that the interactions between nucleons have 
significant noncentral parts. This has several consequences. The magnetic moment of a quasiparticle is not equal to 
the free space value (as discussed above, the same holds for the axial coupling) and it is a tensor, that depends on 
the orientation of the spin with respect to the momentum of t he q uasiparticle. In addition, the Landau quasiparticle 
interaction contains tensor and other noncentral contributions [23j , which couple spin and orbital degrees of freedom. 
For the response to a magnetic field, which is completely equivalent to the case of an axial-vector probe, these effects 
have been explored in Ref. (24[. 

In Landau Fermi-liquid theory, one describes the long-wavelength, low-frequency response of the system in terms 
of quasiparticles. However, if the current of interest is not a conserved quantity, the corresponding response function 
at long wavelengths contains contributions that cannot be expressed in terms of quasiparticle degrees of freedom. 



In addition, there are two-body contributions to the effective operators. In Ref. 11(, it was shown from sum-rule 
arguments that the contribution to the response not coming from single particle-hole pairs could be substantial. One 
class of processes that can be calculated within Landau Fermi-liquid theory corresponds to creating a single particle- 
hole pair, which subsequently creates a second pair. This is taken into account by including a collision term in the 
transport equation for quasiparticles, and in Ref. 12] it was described how to do this, starting from diagrammatic 
perturbation theory. 
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The general formalism for calculating the rates of kinetic processes from microscopic theory is well developed, but 
to apply it to specific physical situations is usually complicated. However, if collisions are sufficiently infrequent, one 
can adopt an approach based on a kinetic equation similar to the Boltzmann equation for dilute gases, in which one 
introduces a distribution function for the elementary excitations that depends on the momentum of the excitation. 
More generally, when the width of an excitation becomes comparable to the real part of the energy of an excitation, 
it is necessary to work in terms of the spectral density for adding a single particle to the system (the imaginary part 
of the single-particle propagator), which is a function of energy as well as of momentum 25, 26]. In this paper, we 
assume that the widths are sufficiently small that a kinetic equation approach can be used. 

Next we describe the quasiparticle transport equation for a single-component Fermi system with spin 1/2. We as- 
sume that the system is not magnetically polarized. The generalization to isospin is straightforward. The quasiparticle 
distribution function is a matrix in spin space and we write it as 

[^p]aa' & aa ' ~\~ Sp • CT aa i . (9) 

Likewise, the quasiparticle energy can be written in the form 

[^plaa' £p &aa' hp • CT aa i , (10) 

where e p and h p are the spin-independent and spin-dependent contributions to the quasiparticle energy. The linearized 
transport equation in momentum space for the spin response 5s p of quasiparticles with momentum p is given by [l^,[22j 

(uj - £p +q /2 + ep-q/2) 8s p + (n p+q/2 - fip-q/2) £h p =il a [s p >] , (11) 

where the perturbation to the quasiparticle energy is 

Sh p = U a + 2 J ^ g pp , Ss p , , (12) 

and the dependence of <5s p (cj,q) and <5h p (oj,q) on the energy and momentum transfers is implicit. Here, /^[sp/] is 
the collision integral, the prime on the momentum argument indicating that it generally depends on the distribution 
function for states other than p, and U CT is an external field that couples to the nucleon spin. The spin-dependent 
Landau quasiparticle interaction has a central part, cri • er 2 , as well as symmetric tensor and antisymmetric 
terms [23j. Since the latter are generally weaker [24 |. we keep only the central term in Eq. (fl2"|) . For the density 
response, Eq. (UJ) holds with the spin-dependent contributions replaced by their spin-independent counterparts, and 
the equation analogous to Eq. (fT2"]) is 

5e p = U + 2j^f pp ,Sn p ,. (13) 

In local equilibrium, the net collision rate vanishes and the distribution function is given by the equilibrium Fermi 
function for quasiparticle energy e p , evaluated at the values of the local chemical potential, temperature, and flow 
velocity corresponding to the local number, energy, and momentum densities. The quasiparticle energy that occurs in 
the local-equilibrium distribution function includes contributions from quasiparticle interactions, so the quasiparticle 
energy is not the one for the equilibrium state. This choice is physically the most meaningful, because in the 
energy conservation condition the quasiparticle energies that appear must include the effect of the non-equilibrium 
quasiparticle distribution. From linear response theory and for u — 0, the local-equilibrium response then follows 
from Eq. (|11[) and is given by 

8s p \ lc = R p 5h p with Rp = n P+q/ 2 Z n P~q/ 2 t (14) 

£p+q/2 - £ p-q/2 

where the subscript "le" denotes the value of the quantity for local equilibrium. 



IV. RELAXATION TIME APPROXIMATION 



In general it is difficult to solve the transport equation for the full collision integral. We therefore approximate the 
collision integral as 

<5s D — Ss D \, 

I a [ Sp ,} = ? ^, (15) 
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where r a is an average relaxation time. In this section, we focus on the spin response, but analogous expressions 
hold for the density and isospin responses. Equation (|15p amounts to the assumption that all angular harmonics of 
the spin-dependent part of the quasiparticle distribution function relax at the same rate, and this form ensures that 
the collision term vanishes when <5s p = <5s p |i c . In addition, the relaxation time is assumed to be independent of the 
quasiparticle momentum. However, consideration of the scattering process in detail shows that, in order to obtain 
agreement with rates in the collisionless limit, |cj|t ct — > oo, the relaxation time must depend on the energy transfer 
(see Ref. [l2| and Sect.fV]). For the spin response, t ct corresponds to the rate of change of the nucleon spin through 
collisions with other nucleons, and by solving the transport equation, we include multiple-scattering effects. 

More generally, one could have allowed for changes in the temperature of the two different spin components, but for 
Fermi systems at low temperatures, this effect, which corresponds to thermoelectric phenomena for charged systems, 
is relatively unimportant. For most condensed matter systems, Eq. (|15[) is a rather poor approximation, since the 
total spin, which corresponds to the component of the deviation function having angular symmetry corresponding 
to I = is conserved to a good approximation because noncentral forces generally play little role, while higher-? 
components of the spin deviation function can decay on a much shorter timescale. For example, in liquid 3 He, the 
lack of spin conservation is due to the interaction between the nuclear magnetic dipole moments, which is very weak 
compared with the central parts of the interatomic interaction. However, in nuclear systems noncentral contributions 
to nuclear interactions, especially those from tensor forces due to pion exchanges, are strong and the single relaxation 
time approximation is expected to be better. The approximate form for the collision term in the transport equation 
for the density response must have a more general form, since particle number conservation ensures that the I = 
component of the distribution function does not relax and, for a single-component system, momentum conservation 
ensures that the I = 1 component does not relax either (see for example Ref. (27j). For a multi-component system, 
such as a mixture of neutrons and protons, the number of particles of each component is conserved, and consequently 
the I = components cannot relax, but the / = 1 components can relax, because momentum may be transferred from 
one component to another. 



A. Calculation of the response function 



With the approximation Eq. (|15[) . the linearized transport equation can be rewritten in the following form 



- v p • q Us p 



R p Sh p = . 



(16) 



with £p+ q /2 — £p-q/2 ~ v p ' q- in the expansion of the quasiparticle interaction in Legendre polynomials, the / = 
term go is the dominant spin-dependent contribution in neutron matter (30| . and therefore we neglect the higher-/ 
terms. With this assumption, the perturbation to the quasiparticle energy, Eq. (|12p . is given by 



Shp = U a + g s with s = 2 



dp' 

(2tt) 3 p ' 



(17) 



As in Eq. ([5]) , s is the Fourier transform of the spin deviation. We then solve the transport equation and find 

s= - X ^,q)U CT , (18) 

where the response function Xa- is given by 



Xo 



1 + g X a 



and X„ 



dp' v p q %jr a 



R, 



(27r) 3 lu + ijr a - v p • q p 



(19) 



Here X a is the response function in the absence of mean-field effects. Provided the temperature is low compared 
with the Fermi energy, the main contributions to the integral in Eq. (|19[) come from the vicinity of the Fermi surface, 
which leads to 



X„ 



N(0) 



1 



cu + i/T a + v F q 
2vpq \bj + ijr n — v-pq 



■In 



(20) 



where N(Q) = m*kp/Tr 2 is the density of states at the Fermi surface for both spin populations, m* being the nucleon 
effective mass and vp — k-p /m* the Fermi velocity. For the imaginary part of \„ we have 



(21) 
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with dimensionless Landau parameter Go = N(0) go, and X a — X a /N(0), whose imaginary part is 



lmX„ = 



2vpq 



arctan[(w + Vyq)^] — arctan[(w — VFq)r a ] 



(22) 



The branch of the arctangent to be used is that lying between —tt/2 and +n/2. For r a — > oo, the form for Im^ given 
by Eqs. ([21^ and (|22|) reproduces the results of Ref. [I[ for single particle-hole pair states, with 



ImX, 



2vpq 



<d(v F q 



>\) 



(23) 



where Q(x) is the step function. Our results generalize earlier work by taking into account effects due to non-zero 
wavelengths and recoil of the nucleons. A direct inspection shows that the resulting dynamical structure factor 
satisfies the detailed balance condition S(— u>) = S(uj)e~ w / T . In contrast to Ref. [12|], where calculations were made 
to leading order in the scattering rate, Eq. (|21l) contains contributions of higher order and thereby takes into account 
the Landau-Pomeranchuk-Migdal effect [28|, |29( . 
In the long- wavelength limit, q — > 0, we have 



X a {u,q^0) 



1 



1 — lUlTcr 



and x CT (o;,g — > 0) 



f 



1 + G - iuJT a 



with imaginary part 



Imx CT (w,<7 -> 0) 



(1 + Go) 2 + K) 5 



(24) 
(25) 



In the absence of mean-field effects, this has the same form as the Ansatz used by Raffelt et al. to account for multiple 
scattering at low ui [Ed, 03. Equation (f2"5| shows that the characteristic frequency for the response is ~ (1 + Go)/r a . 
The factor 1 + Go indicates that near the transition to a ferromagnetic state, Go — > — 1, the characteristic time 
becomes long, corresponding to what is referred to as critical slowing down. For neutrons, one has Go > (30j and 
the spin response is pushed to higher frequencies. 



V. RELAXATION TIMES 



To begin, we consider the time for an excess population of quasiparticles in a particular momentum, energy and spin 
state (denoted by pi, e% and <Ti) to relax when the distribution function for all other states is that for equilibrium. 
It is convenient to consider the general case when the quasiparticles of the excess population are not on the en ergy 
shell, since this is the quantity which naturally enters calculations of the response functions at high frequency [12j . 
The relaxation time can be written in operator form 



1 



f 



t(ei + U), <T\ ■ pi) r(£i+Cj) 



{1 + acTf pi) , 



(26) 



where a is a coefficient that characterizes the strength of noncentral contributions to the relaxation rate. Unlike in 
systems with only central interactions (a = 0), the relaxation rate depends on the spin orientation of the quasiparticle, 
because spin and momentum are coupled. 



By generalizing the standard theory of relaxation rates 22j to the case of noncentral interactions, we have [1 



r(ei 



C 



(ei +w) 



2 1 



where the factor 3/4 is included so that energy-averaged relaxation rates have a simple form (see Eqs 
and the coefficient G is given by 

47T 3 -I— r / TO* 



and 



(27) 

m 



c = 



3JV(0) 2 



^3 S( Pi -k F )) (2ir) 3 6( Pl +p 2 -P3-P4) \ Tr [ , CT2 (k, k') A„ t , CT2 (-k, k') ] 



(27T) 



(28) 



Here we have taken pi to lie on the Fermi surface, yt CTl!(T2 (k, k') denotes the quasiparticle scattering amplitude in 
units of the density of states, k = pi — P3 and k' = pi — P4 are the momentum transfers, 2 and we have neglected 



2 We use k and k' for the momentum transfers between nucleons, in order to distinguish them from the momentum transfer q in the 
structure factors. This differs from the notation used in Refs. 23, 30, 31] and these should also not be confused with relative momenta. 
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the neutrino momenta in the delta function that expresses momentum conservation, since they are small compared 
with the Fermi momentum. The factor 1/4 in Eq. (|28[) is the symmetry factor. 3 Since we work with antisymmetrized 
amplitudes one factor of 1/2 is necessary to avoid double counting of final states, and a second factor of 1/2 comes 
from taking the average over initial spin states of particle 1. On the Fermi surface, the momentum transfers are 
orthogonal and we can express Eq. ([28)1 as 



C 



6fc 



— x iTr[^ CTl , CT2 (k,k')^ CTl , CT2 (-k,k')] ), 



(29) 



where the average is over the Fermi surface. In terms of fc, k' , this can be written as [3ll | 

'dk' fc F e(4fc| -fc 2 -k' 



[F(k,k') 



2kp 2kp 

1 f dk r u ~ <^<ai- 2 — i- 2 — i~' 2 ^ 



7T J fcp 




k F ^/Ak 2 -k 2 -k'< 



F(k,k') 



(30) 



With this average, the coefficient a can be written as 



a= 2 



1 / Tr [ o-i • pi A ai ltT2 (k, k') A ai , CT2 (-k, k') ] 



Tr[^l CTl:CT2 (k,k')^ CTl , CT2 (-k,k')] 



(31) 



More general disturbances of the quasiparticle distribution function will depend both on the direction of the quasi- 
particle momentum on the Fermi surface and on the spin of the quasiparticle, and the relaxation time for the distur- 
bance will depend on an average of the scattering rate over the Fermi surface and over quasiparticle spins, weighted 
by functions of the direction of the quasiparticle momentum and of the spin. In general, the eigenstates of the col- 
lision operator will have a definite value of the total angular momentum, which is made up of an orbital component 
coming from the dependence of the quasiparticle distribution on the angle on the Fermi surface and of the spin of the 
quasiparticle. 

The most important case for relaxation of long-wavclcngth spin fluctuations is a disturbance of the distribution 
function corresponding to a spin polarization that is independent of direction on the Fermi surface. For long wave- 
lengths \u>\ ^> vpq and for frequencies large compared with the relaxation rate \u>\ ^ 1/tv, the appropriate average 
relaxation time for the transport equation and the spin response is given by [12j 



1 



1 



ujN(0) 



E 

m s1 



dpi n{e\) — n{e\ + uj) 



(27r) 3 T a (ei + W,<Ti • pi) 



(32) 



^ 1 f n(ei) - n(ei +u>) 

UJ J 1 T a (e!+U}) 



(33) 



where the noncentral term in the spin relaxation rate (a CT in the operator form analogous to Eq. (|26[) ) averages to 
zero. Following Refs. [l2|, HH, one has for the coefficient C a for the spin relaxation rate 



6fc| 



J=l,2,3 



Ar 1 ,«r 2 (k,k')o^[(<Ti +tT2 y t A rTl , CT2 (-k, k')] 



(34) 



The commutator with the two-body spin operator demonstrates that only noncentral terms in the scattering amplitude 
contribute. The factor 1/12 in Eq. (|34| includes the symmetry factor 1/4 and a factor 1/3, because we have summed 
over all possible directions of the spin component j. 

Since the dependence on the quasiparticle energy factorizes from the nuclear interaction part, we can directly 
calculate the average relaxation time of Eq. (|33[) and finally obtain 



C [T 2 + (u/2ir) 2 



and 



C a [T 2 + (lu/2it) 2 



(35) 



3 We note that Refs. [I3I, |3ll use a symmetry factor of 1/2 instead of 1/4 and consequently overestimate rates by a factor 2. 
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A. One-pion exchange interaction 



For the one-pion exchange (OPE) interaction, the direct and exchange contributions to the scattering amplitude in 
Born approximation are given by 



<^ 2 (k,k')=-iV(Q) 



9a 

2F« 



2 r 



o"! ■ ker 2 • k <Ti • k' ct 2 • k' + k' 2 (l - <7, ■ er 2 )/2 



k 2 + ml 
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(36) 



with pion decay constant = 92.4 MeV and neutral pion mass — 134.98 MeV. The spin trace in the relaxation 
time for the spin response, Eq. (|3~4")1 , leads to 
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(37) 



For ttt-tt = 0, each of the three terms in the square bracket of Eq. (|3T[) yields 1 when averaged over the Fermi surface 
according to Eq. (|30p. and for non-zero m w this integral can be calculated analytically, and one finds for the spin 
relaxation rate from one-pion exchange fl3j 
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OPE 
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where the factor G{x) takes into account the effects of a non-zero pion mass 
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(38) 



(39) 



For \uj\t,, 3> 1, the imaginary part of the spin response function in the long- wavelength limit is given by N (0) / ((VT a ) 
(see Eq (|25|) ). In this limit, when multiple-scattering effects are small, our result for the dynamical structure factor 
using the spin relaxation time of Eq. (|3"8|) agrees with the result of Raffelt et al. [!,[§] using f/m^ fa g a j2F^. 

We can compare the spin relaxation time r° PE with the relaxation time corresponding to decay of an excess of 
quasiparticles in a particular momentum state t . For the latter, the spin trace of Eq. ([37]) has to be replaced by 
the one in the brackets (. . .) of Eq. (|2"9")l. which yields exactly the same result as the right-hand side of Eq. (|3"T|) up 
to the factor 4/3. As a result, we find that the spin relaxation rate and thus spin- flipping collisions obtained from 
the one-pion exchange interaction are comparable to the relaxation rate for decay of an excess population in one 
momentum state, with 



r OPE 



3 T 



OPE 



(40) 



This result highlights the importance of noncentral contributions to nuclear interactions and encourages us to perform 
more systematic calculations of these rates beyond one-pion exchange. Next, we calculate the contributions to the 
relaxation times from a general representation of the quasiparticle scattering amplitude and present results in Sect. I VII 



B. General operator representation 



For neutron matter, using the general operator representation of the scattering amplitude on the Fermi surface in 
the notation of Refs. [23, bLU], we find for the spin trace of Eq. (f34|) : 4 



- Y Tr 
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. (41) 



4 We note that the factor 3 in front of the cross vector amplitude in Eq. (7) of Ref. [3lll should be 1. 
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where the amplitudes on the right-hand side are functions of k and k' . The scattering amplitudes on the Fermi sur- 
face Aonsor, -Aexch. tensor, Apin-orbit , Aiiff. vector and .Across vector are real and characterize the momentum-dependent 
strengths (in units of the density of states) of the tensor operator 5i2(k), the exchange tensor Si2(k'), the spin-orbit 
operator i{cr\ + <r 2 ) • k x k', the spin difference vector i(<xi — <r 2 ) • k x P (or antisymmetric spin-orbit), and the cross 
vector operator (<Ti x go) ■ (k' x P), respectively, with two-body center-of-mass momentum P = Pi + P2 = Pa + P4 
(for details, see Refs. [23jlSl|)- The latter two operators do not conserve the spin of the interacting particle pair and 
are induced in the medium due to screening by particle-hole excitations [23[ . Finally, the tilde on the tensor parts of 
the scattering amplitude indicates that they take into account induced center-of-mass tenso r op erator contributions, 
since this is not a linearly- independent operator on the Fermi surface, as discussed in Refs. [23l. l3l|. 

For the spin trace of Eq. (|29[) corresponding to the relaxation rate for decay of an excess population in one momentum 
state, we have 



Tr [ Am ,ct 2 Am ,ct 2 ] A 
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(42) 



where in addition the central parts of the scattering amplitude, *4 SC aiar and -4 sp i n , contribute. These correspond to 
the spin-independent amplitude and the spin-spin operator cri • cr 2 , respectively. We note that all contributions in 
Eqs. (HH) and (|42p are positive. The minus sign of the direct-exchange tensor interference term is canceled by a 
relative minus sign in the exchange tensor amplitude. 



VI. RESULTS 

We calculate the contributions beyond one-pion exchange based on low- momentum interactions V\ ow k [H|, |33| |. 
which are obtained by evolving nuclear forces to low momentum using the renormalization group. The resulting two- 
nucleon interactions become universal at momentum scales A < 2fm _1 and provide a basis for model-independent 
predictions of low-energy processes. The renormalization-group evolution preserves the long-range parts from pion 
exchanges and V\ ow k includes subleading noncentral contributions, so that all low-energy nucleon-nucleon scattering 
observables and deuteron properties are reproduced. In this first study, we have not included contributions from 
low-momentum three- nucleon interactions [341 ] . Their effects are generally weaker in neutron matter, but calculations 
of the equation of state show that three- nucleon interactions become important for fcp ^ 1.5 fm -1 [351 ]. We will study 
their contributions to neutrino processes in future work. 

In addition, we include many-body noncentral and central correlations from second-order particle-particle (plus 
hole-hole) and particle-hole contributions using the same Viowfc interactions. The resulting quasiparticle scattering 
amplitudes are discussed in detail in Ref. |23l| and have been used to calculate the neutrino emissivity from pair 
bremsstrahlung for neutron star cooling 31]. Based on our results and general arguments [36| . second-order cor- 
rections become reasonable for low-momentum interactions. The intermediate states include all possible excitations 
for interacting particles on the Fermi surface. We use the effective mass obtained from the lowest-order Vi ow fe for 
all results, including for the estimates based on the one-pion exchange interaction. The effective mass varies from 
m*/m = 0.95 at &f = l.Ofm -1 to m*/m — 0.78 at fcp = 2.0fm~ 1 , and in this range is well approximated by a linear 
dependence on the Fermi momentum. We note that one expects an increase of the effective mass due to polarization 
effects, but this is compensated by the reduction of the quasiparticle strength Zfc as can be seen from the results of the 
renormalization-group calculation of induced interactions in neutron matter [30| ■ We emphasize that a second-order 
calculation cannot give final results, but it provides a range for the effects due to many-body correlations. 

Finally, we note that the effect of particle-particle correlations on neutrino-pair bremsstrahlung and other neutrino 
processes has been investigated previously in Refs. [37l l38l l39l. EoT] . 



A. Relaxation times 



Our results for the spin relaxation coefficient C a of Eq. (|34[) are shown in Fig. Q] For energies u> = and T ~ 
5 — lOMcV, the value of C a = O.lMeV -1 corresponds to spin relaxation rates \/r a — 2.5 — lOMeV. We find that 
the OPE model significantly overestimates the strength of noncentral contributions, compared to low- momentum 
interactions Vlowfc, for all considered densities. Beyond the V\ ow k results, we find that second-order many-body 
contributions reduce the spin relaxation rate especially at lower densities (note that C a is proportional to the square 
of the quasiparticle scattering amplitude) . These effects are due to second-order particle- hole interference of tensor 
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FIG. 1: (Color online) The spin relaxation rate given by C CT of Eq. (|34|l as a function of Fermi momentum &f obtained from 
the one-pion exchange interaction (OPE), from low-momentum interactions Viowfc, and including second-order many-body 
contributions. In addition, we show that the result obtained from VW k plus second-order contributions is dominated by tensor 
interactions (dotted versus solid line). 




FIG. 2: (Color online) The relaxation rate for decay of an excess of quasiparticlcs in a particular momentum state given by C 
of Eq. (|29|l as a function of Fermi momentum kp obtained from the one-pion exchange interaction (OPE), from low- momentum 
interactions Viowfc, and including second-order many-body contributions. In addition, we show that the result obtained from 
Viawh plus second-order contributions is dominated by central interactions (dotted versus solid line). 



with strong central interactions, which are driven by large scattering lengths at very low densities. The band in Fig.[T] 
from Viowfc to including second-order contributions provides a range for the effects due to many-body correlations. In 
addition, we observe that the spin relaxation rate depends only weakly on density, and the rate obtained from V\ ow f. 
plus second-order contributions is dominated by the tensor terms in Eq. (|4ip . 

For the relaxation coefficient C of Eq. corresponding to decay of an excess of quasiparticles in a particular 
momentum state, we obtain rates in Fig. [5] that are of similar magnitude compared with the spin relaxation rate. 
While the OPE rate is approximately independent of density, the OPE model underestimates the relaxation rate at low 
densities. This is because the central part of the OPE interaction ~ k 2 and ~ k' 2 does not capture the central shorter- 
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FIG. 3: (Color online) Ratio of the spin relaxation rate to the relaxation rate for an excess of quasiparticles in a single 
momentum state (1/t ct )/(1/t) as a function of Fermi momentum &f for purely tensor scattering amplitudes (in which case 
the value is 2), for the one-pion exchange interaction (which gives the value 4/3), from low-momentum interactions V\ ov ,k, and 
including second-order many-body contributions. 




CO/(v p q) 



FIG. 4: The imaginary part of the spin response function lmx^/N(0) of Eq. (|21[1 in units of the density of states versus u}/(vpq). 
Results are shown for the non-interacting system, without and with mean-field effects, Go = and Go = 0.8 respectively, and 
for different values of the spin relaxation rate l/r CT = 0, VFqr a = 2 and VYqr a = 5. 



range physics in nuclear forces. This deficiency of the OPE model is most prominent at low densities, in comparison 
to the increasing Vj ow k rate. Similar to the spin response, we find a reduction of C due to second-order many-body 
contributions, where the band in Fig. [5] again indicates a range for the effects due to many-body correlations. Finally, 
as expected, the relaxation rate obtained from V\ ow & plus second-order contributions is now dominated by the central 
terms in Eq. (j42j) . 

In Fig. [3] we show the ratio (l/r cr )/(l/r) of the spin relaxation rate to the relaxation rate for an excess of quasi- 
particles in a single momentum state as a function of Fermi momentum fcp- This is a very useful measure of the 
strength of noncentral interactions compared to central ones. For purely tensor scattering amplitudes, the ratio of 
the corresponding spin traces in Eqs. (|4"Tj) and (|42j) gives (1/t ct )/(1/t) = 2, while for the OPE interaction, which 
has a central part in Eq. ([36]) . this ratio is (I/t^/^/t) =4/3, see Eq. (|40|) . While the ratio obtained from Mowfc 
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Go 


0.8 


0.8 


0.8 


k F [fin" 1 ] T [MeV] 


C a from OPE 


Viow k 


Mow k + 2nd order 


1.0 5 
10 


0.0770 0.0697 
1.08 0.798 


0.0397 0.0386 
0.612 0.554 


0.00754 0.00753 
0.120 0.120 


1.7 5 
10 


0.119 0.107 
1.66 1.21 


0.0476 0.0468 
0.744 0.700 


0.0296 0.0294 
0.470 0.457 



TABLE I: Thermally averaged inverse neutrino mean free path (A -1 ) in km -1 calculated from Eq. (|44|l for characteristic 
temperatures and Fermi momenta. Results are given without and with mean-field effects, Go — and Go = 0.8 respectively, 
and for different spin relaxation rates l/r CT based on Fig. [1] 

and including second-order many-body contributions is considerably smaller at low densities, the relative strength of 
noncentral interactions increases with momentum and thus with density, as can be seen in the results of Fig. [3] based 
on modern nuclear forces. 

B. Dynamical structure factor 

Motivated by the importance for neutrino rates, we focus on the spin response in this section. The dynamical 
structure factor is determined by the imaginary part of the spin response function Imx <x, which is given by Eq. (|2ip in 
the relaxation time approximation. In units of the density of states, the imaginary part ImXa is a function of vpqT a and 
to/ (vpq) or of vpqT a and u>r a - In the long-wavelength limit, q — > 0, we have already found that this is proportional to 
uj times a Lorentzian function of to, see Eq. (|25[) . Therefore, we plot in Fig. Uthe imaginary part of the spin response 
function versus u)/(vpq). Results are shown for the non-interacting system, without and with mean- field effects, 
Go = and Go = 0.8 respectively, and for different values of the spin relaxation rate l/r CT = 0,vpq/5, and vpq/2. 
We have taken the Landau parameter from renormalization-group calculations of induced interactions in neutron 
matter [3(3], which yield Go ~ 0.8 over the densities considered in Sect. IVI Al The values of Vpqr a = 2 — 5 correspond 
to spin relaxation rates based on Fig. [T]for typical momentum transfers q ~ lu over the range T = 5 — 10 MeV and 
kp = 1.0 — 1.7fm _ . With l/r CT comparable to vpq, these estimates also show that recoil effects may be important. 

In the non-interacting case, Go = and l/r CT = 0, the imaginary part of the spin response function is given by 
■KLo/(2vpq) times a step function, see Eq. (f2"3"| . With single-pair mean-field effects, Go = 0.8, a collective spin-zero- 
sound mode appears as a pole contribution at u> / '(vpq)\ zs > 1, where the position of the pole is given by Q 

l + GoX a (u/(v F q)\ zs ,l/T a = 0) = 0. (43) 

As the spin relaxation rate increases, going from l/r CT = to vpq/5 and vpq/2 : the response is pushed to higher 
frequencies and the spin-zero-sound peak disappears already for these moderate spin relaxation rates. For comparison, 
we also show the effects due to single-pair states at vpqr a = 5, where interactions (Go = 0.8) decrease the response 
at low u)/(vpq) and also move the strength to higher frequencies. 

C. Neutrino mean free paths, energy loss and energy transfer 

We next assess the significance of the improved rates for neutrino mean free paths, energy loss and energy transfer. 
For derivations of Eqs. (|4"4")) to (|4"T)) see Refs. All rates are for one neutrino flavor. We emphasize that the 

OPE results are based on the solution to the transport equation in the relaxation time approximation, and do not 
correspond directly to OPE rates used in supernova simulations. For simple estimates, we use the dynamical structure 
factor for spin fluctuations in the long- wavelength limit, S a {uj) = S a (ui,q — > 0), given by Eqs. and (|2"B")) . without 
further approximations or Ansaetze for the structure factor. Effects due to the finite wavelength and recoil of the 
nucleons will be studied in future work. 

In Table U we present results for an average inverse neutrino mean free path (A -1 ), 

oo 

(rl> = ^J/ 4,w,e " /TWI (44) 



for characteristic temperatures and Fermi momenta. This result applies for a Maxwellian initial distribution of 
neutrinos, and Pauli blocking in the final state has been ignored. We consider structure factors without and with 
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Go 


0.8 


0.8 


0.8 


k F [fin" 1 ] T [MeV] 


C a from OPE 


Mow k 


Mow k + 2nd order 


1.0 5 
10 


1.77 1.62 
4.02 3.00 


0.911 0.888 
2.25 2.06 


0.173 0.172 
0.441 0.440 


1.7 5 
10 


2.75 2.49 
6.18 4.55 


1.09 1.07 
2.73 2.57 


0.679 0.675 
1.72 1.68 



TABLE II: Energy-loss rate Q of Eq. (|45[) due to neutrino-pair bremsstrahlung, nn — > nnvv, for characteristic temperatures 
and Fermi momenta. Results are given without and with mean-field effects, Go = and Go = 0.8 respectively, and for 
different spin relaxation rates 1/t ct based on Fig. [1] The energy-loss rates are in units of 10 33 erg cm -3 s" 1 for T = 5 MeV and 
10 35 ergcm _3 s- 1 for T = 10 MeV. 



Go 





0.8 





0.8 





0.8 




fc F [fm- 1 ] 


T [MeV] 


C a from OPE 


How k 


VJowfe - 


- 2nd order 






5 


2.48 


2.26 


1.27 


1.24 


0.241 


0.241 


nn *- 


> nnvV 


1.0 


3.46 


2.81 


1.94 


1.76 


0.401 


0.394 


vnn 


—> vnn 




10 


2.81 


2.10 


1.58 


1.44 


0.308 


0.307 


nn *- 


> nnvv 




3.41 


2.24 


2.20 


1.79 


0.502 


0.485 


vnn 


— > vnn 




5 


3.85 


3.48 


1.53 


1.50 


0.949 


0.943 


nn *- 


> nnvv 


1.7 


5.33 


4.30 


2.38 


2.20 


1.53 


1.46 


vnn 


— > vnn 




10 


4.32 


3.18 


1.91 


1.80 


1.21 


1.18 


nn *- 


> nnvv 




5.21 


3.37 


2.76 


2.35 


1.84 


1.67 


vnn 


—> vnn 



TABLE III: Rate of energy transfer AQ/AT due to neutrino-pair bremsstrahlung and absorption, nn <-> nnvv, of Eq. (|46[) 
and due to inelastic scattering, vnn <-> vnn, of Eq. (|470 for characteristic temperatures and Fermi momenta. Results are given 
without and with mean-field effects, Go = and Go = 0.8 respectively, and for different spin relaxation rates l/r a based on 
Fig. [1] The rates are in units of 10 33 erg cm" 3 s _1 MeV -1 for T = 5 MeV and 10 35 erg cm" 3 s" 1 MeV -1 for T = 10 MeV. 

mean-field effects, Go — and Go = 0.8 respectively, and for different spin relaxation rates l/r a based on Fig. [T] 
With the spin relaxation rates obtained from Vi ow k and including second-order many-body contributions, the mean 
free paths are significantly longer compared to the OPE model. This follows the reduction of C a seen in Fig. [TJ For 
OPE, the effects of interactions (Go = 0.8 compared to Go = 0) reduce the neutrino scattering rate, especially at 
higher temperature. In contrast, with the rates based on low-momentum interactions, Lor a is larger and the imaginary 
part of the spin response function approaches Imx CT (o;, q — > 0) — * N (0) / (uiT a ) . As a result, mean-field effects arc weak 
for \oj\T a ^> 1 in the long-wavelength limit. 

The energy-loss rate Q due to neutrino-pair bremsstrahlung, nn — > nnvv, of neutron matter transparent to neutrinos 
is given by 

oo 

Q= G|G|n r^e-u/Tg^y (45) 

MIT J 



Our results for the energy- loss rate Q arc listed in Table [TT] for characteristic temperatures, Fermi momenta and 
the different cases of the structure factor. They follow the same general pattern as the inverse mean free paths in 
Table HI a reduction of the energy loss calculated with modern nuclear forces compared to OPE and consequently 
weak mean-field effects. 

Finally, we consider the rate of energy transfer AQ/AT from neutron matter at temperature T to a neutrino fluid 
at temperature T„, with AT = T — T u and |AT| <C T. The energy transfer due to neutrino-pair bremsstrahlung and 
absorption, nn <-* nnvv, is given by 

o 
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and for inelastic scattering, vnn <-> i/nn, one has 



OO 



AQ _ 30C 2 A G 2 nT 3 
AT ~ 10^3 



/ 



duj ui 2 (12 + Quj/T + (uj/T) 2 ) e 



(47) 



o 



Our rates for the energy transfer are shown in Table IIIII for characteristic temperatures, Fermi momenta and the 
various cases for the structure factor. The pattern of these rates is similar to what we found for other rates in Tables U 
and [Til In addition, for all cases we find that the energy transfer due to inelastic scattering is less than a factor 
2 larger compared to the contributions from neutrino-pair bremsstrahlung and absorption. In contrast, Hannestad 
and Raffelt estimated this ratio to be 10 [10] . However, in making this estimate they used Eqs. ([4"1?[) and (|4"T[) with 
IiriXo-(uj,q — > 0) ~ 1/w 2 , while we find Imx CT (w,q — > 0) ~ l/u>. 



We have developed a unified treatment for neutrino processes in nucleon matter based on Landau's theory of Fermi 
liquids that includes one- and two-particle-hole pair states consistently. The contributions from two-particle-hole pair 
states are crucial for neutrino-pair bremsstrahlung and absorption, for inelastic scattering, modified Urea reactions, 
and axion emission. In supernovae, neutrino-pair bremsstrahlung and absorption dominate the neutrino-number 
changing reactions and are key to the production of muon and tau neutrinos 

Neutrino rates involving two nucleons can be calculated in terms of the collision integral in the Landau transport 
equation for quasiparticles. Using a relaxation time approximation, we have solved the transport equation for density 
and spin-density fluctuations and derived a general form for the response functions. The solution includes multiple- 
scattering effects and effects due to non-zero wavelengths and recoil of the nucleons. We have applied our approach 
to neutral-current processes in neutron matter, but the generalization to isospin is straightforward. Our results for 
the spin response are summarized by Eqs. ©, (12"0[) , (|2"Tj) . (l2"2"j) . ([33)1 and the values of C a of Fig.[TJ 

We have calculated the relaxation times based on the OPE model and for a general representation of the quasiparticle 
scattering amplitude. For OPE, the spin relaxation rate is comparable to the quasiparticle relaxation rate, r /t„ — 4/3. 
This highlights the importance of noncentral contributions to nuclear interactions. We therefore performed more 
systematic calculations of these rates. In addition, for \u>\r a ^S> 1 and in the long-wavelength limit, our result for the 
dynamical structure factor agrees with Raffelt et al. [1, [§] . 

Beyond OPE, we have calculated the relaxation times based on low-momentum interactions V\ ow k and including 
second-order many-body contributions. The effects of three-nucleon interactions are generally weaker in neutron 
matter [351 ]. but need to be included in future work. The OPE model significantly overestimates the strength of 
noncentral contributions, compared to low-momentum interactions Vj ow fc, for all considered densities. Beyond the 
V\ow k results, we have found that second-order many-body contributions reduce the spin relaxation rate especially at 
lower densities. This provides a range in Figs. [1] and [2] for the effects due to many-body correlations. By using spin 
relaxation times that incorporate both "in-scattering" and "out-scattering" terms in the transport equation, effects 
corresponding to vertex corrections in the microscopic theory are automatically taken into account. 

Using the spin response in the long- wavelength limit, but without further approximations or Ansaetze for the 
structure factor, we have estimated the significance of the improved rates for neutrino mean free paths, energy loss 
and energy transfer. We have found a reduction of these rates using modern nuclear forces compared to OPE and 
consequently weak mean-field effects. In addition, for all cases we find that the energy transfer due to inelastic 
scattering is not significantly larger than that due to neutrino-pair bremsstrahlung and absorption. 

One may ask how good the relaxation time approximation is. Our choice of spin relaxation time is designed to 
agree with microscopic theory in the collisionlcss limit, \<jj\t 1, and at long wavelengths. For the hydrodynamic 
limit, \u)\r a <C 1, and long wavelengths, exact solutions of the transport equation have been obtained, and one 



The relaxation time for the hydrodynamic limit is always greater than or equal to that for the collisionless limit. 
For Tct/t ^> 1, r CT |h y dro = T <r, while for t^/t — 1, 7V | hydro = {it 2 /Q)T a . Since for realistic nuclear interactions, r CT 
is significantly larger than r, this indicates that differences between spin relaxation times in the collisionless and 
hydrodynamic limits are expected to be of order a few per cent. Consequently, uncertainties due to the use of the 
relaxation time approximation are small compared with other uncertainties in the calculation. 



VII. CONCLUDING REMARKS 





(48) 
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The use of the quasiparticle transport equation with a collision term allows us to include some two-particle-hole 
pair states, but not all. Among contributions not included are terms that correspond to the incoherent parts of the 
propagator for a particle-hole pair, that is to contributions that do not correspond to an intermediate state containing a 
well-defined quasiparticle together with a well-defined quasiholc. Moreover, there are intrinsic two-body contributions 
to hadronic weak currents. Further work is needed to determine how important these additional contributions are. 

There are numerous directions for future work. One is to explore mixtures of neutrons and protons. A second 
is to extend the calculations to situations when matter is less degenerate. As one sees from our results, there is 
significant uncertainty in the effects of the medium on quasiparticle scattering amplitudes, since there are sizable 
differences between rates obtained with Vl ow /c and those that include many-body contributions to second order, and 
an important task is to reduce these uncertainties. 
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